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Abstract 


This  report  addresses  the  initial  phase  of  hypersonic  boundary-layer  transition  con^rising 
excitation  of  imstable  normal  modes  (receptivity  problem)  and  their  downstream  evolution 
from  receptivity  regions  to  the  instability  growth  onset  (instability  pre-history  problem).  It  is 
shown  that  this  phase  is  more  complicated  than  in  subsonic  and  moderate  supersonic  cases 
due  to  the  following  features  of  the  disturbance  field:  (1)  the  first  and  second  modes  are 
synchronized  with  acoustic  waves  near  the  body  nose  region;  (2)  fiirther  downstream  the 
first  mode  is  synchronized  with  entropy/vorticity  waves;  (3)  near  the  instability  growth  onset 
the  first  mode  is  synchronized  with  the  second  mode.  Disturbance  behavior  in  the 
synchronism  regions  (2)  and  (3)  are  studied  using  the  multiple-mode  method,  which 
accounts  for  interaction  between  modes  of  discrete  and  continuous  spectrum  due  to 
nonparallel  effects  of  the  mean  flow.  It  is  shown  that  vortidty/entropy  waves  are  partially 
swallowed  by  the  boundary  layer  and  effectively  generate  the  first  mode  due  to  Synchronism 

(2) .  This  mechanism  can  conq)ete  with  the  leading  edge  receptivity  to  the  freestream 
acoustic  waves  in  cases  of  “quiet”  freestream  conditions  and  conical  body  configurations. 
Spectrum  topology  and  interaction  between  the  first  and  second  modes  due  to  Synchronism 

(3)  is  studied  using  asymptotic  and  numerical  analyses.  Branch  points  of  the  discrete 
spectrum  were  identified.  It  is  shown  that  near  the  branch  points  normal  modes  have  singular 
behavior.  The  inter-mode  exchange  rule  coupling  input  and  output  characteristics  of  the  first 
and  second  modes  crossing  the  branch  point  vicinity  was  established.  Combination  of  the 
receptivity  estimates  related  to  Synchronism  (2)  and  the  inter-mode  exchange  rule  related  to 
Synchronism  (3)  allows  the  evaluation  of  instability  initial  anq)litudes  required  for  the  PSE 
calculations  of  the  transition  onset  point.  That  makes  feasible  an  anq)htude  method  coupling 
receptivity  and  stability  problems.  Results  of  analyses  are  conq)ared  with  e?q)erimental 
observations  in  the  AEDC  Tunnel  B  on  sharp  cones  at  the  freestream  Mach  number  M  =  8 . 
Possible  ways  of  experimental  investigation  of  hypersonic  receptivity  and  instability  pre¬ 
history  are  discussed. 
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Nomenclature 
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x  =  x*  /L* 
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mass  flux  disturbance  anq)htude 

pressure  disturbance 

Reynolds  number  based  on  the  boundary-layer  scale  ^|v]x*  / 17* 

Reynolds  munber  based  on  the  boundary-layer  scale  -y/xV*  /t/* 

mean-flow  ten^erature 
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mean-flow  streamwise  and  vertical  velocity 
disturbance  of  x  -  and  y -velocity  con^onent 

streamwise  coordinate 
coordinate  normal  to  the  wall  surface 
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s  =  \IRl 

7]  =  y/yfx 
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V 
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co  =  co*6IUl 
'F 


small  parameter 
boundary-layer  coordinate 
viscosity 

kinematic  viscosity 
temperature  disturbance 
Prandtl  nmnber 
angular  frequency 

disturbance  vector 


Subscripts 
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cr 
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in 
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adiabatic  wall 

synchronism  between  entropy/vorticity  waves  and  the  first  mode 

branch  point 

critical 

upper  boundary  layer  edge  or  entropy  wave 

imaginary  part 

initial 

real  part 

vorticity  wave 

wall 


Superscripts 

bar  complex  conjugate  value 

asterisk  dimensional  value 


Abbreviations 

PSE 

Im 

Re 


parabolized  stability  equations 
imaginary  part 
real  part 
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1.  Introduction 


Laminar-turbulent  transition  in  a  hypersonic  boundary  layer  causes  the  significant  increase  in 
heat  transfer  and  skin  fiiction.  Since  the  aerodynamic  coefficients  of  hypersonic  vehicles  are 
sensitive  to  these  changes,  the  abihty  to  predict  the  transition  location  and  simulate  the 
transition  zone  is  uuportant  in  the  design  and  optimization  of  aerospace  vehicles  [1].  The 
stabihty,  control,  and  structural  design  are  also  affected  due  to  the  increased  thermal  and 
aerodynamic  loading.  Recent  technology  does  not  provide  accurate  prediction  of  transition 
due  to  the  complexity  of  this  phenomenon.  Uncertainty  in  transition  location  leads  to 
diminished  vehicle  performance,  primarily  because  of  the  additional  weight  of  the  thermal 
protection  [2].  These  factors  stimulate  fiuther  theoretical  and  experimental  studies  of 
hypersonic  transition  [3]. 

For  flight  and  quiet  wind  tunnel  conditions,  fieestream  disturbances  and  perturbations 
generated  by  hypersonic  body  surfaces  are  normally  small  in  a  frequency  band  related  to  the 
boundary-layer  transition  process.  In  this  case,  the  initial  phase  of  transition  comprises 
excitation  of  unstable  normal  modes  (receptivity  problem),  their  downstream  amplification 
(stability  problem)  and  non-linear  mechanisms,  which  trigger  transition  at  critical  anq)litudes 
of  the  boundary-layer  disturbances  (transition  onset  problem)  [4],  [5]. 

Linear  stability  theory  allows  the  identification  of  unstable  modes  and  prediction  of  then- 
downstream  growth  [3],  [6],  [7].  Stability  theory  is  well  developed  for  subsonic  and 
moderate  supersonic  boundary  layers  [8].  Its  results  are  consistent  with  many  experiments 
[9]-[ll].  Good  progress  has  been  recently  made  in  receptivity  at  subsonic  speeds  [12]-[18]. 
Combination  of  the  receptivity  theory  with  stabihty  analysis  based  on  the  parabolized 
stabihty  equations  (PSE)  [19],  [20]  makes  feasible  an  amphtude  method  [21]  coupling  the 
external  disturbance  characteristics  with  the  transition  onset  locus,  l^uch  progress  toward 
the  subsonic  and  moderate  supersonic  transition  prediction  is  partially  due  to  the  fact  that 
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the  boundary-layer  disturbance  spectrum  is  relatively  single  in  these  cases,  l.e.  unstable 
normal  modes  of  discrete  spectrum  has  no  branch  points,  turning  points  and  other 
singularities.  For  example,  transition  on  subsonic/supersonic  cones,  flat  plates,  axisymmetric 
bodies  at  zero  angle  of  attack  and  airfoils  is  due  to  excitation  and  dovmstream  amplification 
of  the  ToUmien-Schlichting  waves,  which  belong  to  the  first  normal  mode  according  to 
Mack’s  classification  [22].  Other  modes  are  stable  and  their  eigenvalues  are  essentially 
different  fi’om  the  first  mode.  In  this  case,  interaction  between  ToUmien-Schlichting  waves 
and  other  normal  modes  can  be  neglected  in  the  instability  growth  prediction.  Such  a  single¬ 
mode  approach  is  currently  used  in  the  -method  [7]  and  the  PSE  analysis  [20]. 

In  contrast  to  moderate  supersonic  flows,  the  initial  phase  of  hypersonic  transition  (at  Mach 
numbers  larger  than  4)  reveals  the  foUowing  new  features: 

(i)  Besides  the  first  mode,  the  second  and  higher  modes  coexist  in  the  boundary  layer.  They 
belong  to  the  family  of  trapped  acoustic  normal  waves.  Once  the  second  mode  sets  in,  it 
becomes  the  dominant  instabihty  since  its  growth  rate  tends  to  exceed  that  of  the  first  mode. 
For  insulated  surfaces,  this  occurs  for  Mach  numbers  M  >  4 .  For  cooled  surfaces,  the 
second  mode  can  dominate  at  even  lower  Mach  number.  The  existence  of  the  second  mode 
instability  was  established  by  the  experimaits  of  KendaU  [23],  Demetriades  [24],  and  Stetson 
et  al.  [25],  [26].  The  Mach  8  stability  and  transition  e5q)eriments  [25]  (see  also  [26])  for  the 
boundary  layer  on  a  sharp  cone  indicated  that  the  xmstable  high-fi-equency  second  mode 
plays  a  major  role  in  the  conical  boundary-layer  transition.  These  data  are  consistent  with  the 
second  mode  stability  calculations  [27]-[29].  The  linear  stabihty  theory  is  in  a  qualitative 
agreement  with  the  conical  boundary  layer  ejqperiments.  However  quantitative  discrepancy 
stiU  exists.  Stetson  and  Kimmel  [26]  noted  that  “for  the  linearized  stability  calculations,  there 
is  concern  about  the  effects  of  the  large  second  mode  disturbances  and  the  early  presence  of 
the  nonlinear  disturbances.  For  the  experimental  data  there  is  concern  about  the  effects  of 
the  uncontroUed  fi-eestream  environment  and  the  limitations  of  the  hot-wire  anemometry 
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data.”  In  addition,  this  discrepancy  may  relate  to  the  disturbance  spectrum  singularities  near 
the  branch  points  discussed  in  Section  4. 

(ii)  The  second  mode  instabihty  is  associated  with  synchronization  between  the  first  and 
second  modes  [30],  [31],  Near  the  synchronism  point,  the  eigenvalue  spectrum  splits  into 
two  branches.  In  the  vicinity  of  branch  points,  the  eigenmode  decomposition  is  not  valid  and 
should  be  replaced  by  a  local  solution.  Asymptotic  analysis  of  this  solution  reveals  a  strong 
exchange  between  modes  due  to  nonparaUel  effect  [32].  This  finding  leads  to  the  conclusion 
that,  instead  of  a  single-mode  stability  analysis  recently  used  in  the  transition  prediction 
technology,  a  multiple-mode  analysis  is  required  to  account  for  interaction  between  normal 
modes,  particularly  between  the  first  and  second  modes. 

(ii)  Stability  and  transition  experiments  on  a  planar  boimdary  layer  [33]  showed  that  the  low 
frequency  disturbances  are  growing  despite  the  fact  that  the  linear  stability  numerical  results 
indicated  the  first  mode  disturbances  should  be  stable  [29].  Stetson  et  al.  [33]  reported: 
“There  was  no  evidence  of  any  second  mode  harmonics  in  the  planar  boundary  layer.  The 
major  disturbances  are  the  low  frequency  disturbances  which  are  growing  in  a  frequency 
band  expected  to  be  stable.”  These  results  are  consistent  with  the  flat  plate  results  of  Kendall 
[23]  obtained  at  Mach  numbers  3.,  4.5,  and  5.6.  Kendall  reported  that  “fluctuations  of  all 
frequencies  were  observed  to  grow  monotonically  larger  in  the  region  of  a  boundary  layer 
extending  from  the  flat  plate  leading  edge  to  the  predicted  location  of  instabihty;  i.e.  in  a 
region  where  no  growth  was  e?q)ected”.  Mack  [34]  developed  a  forcing  theory,  which  was 
successfully  apphed  to  Kendall’s  Mach  4.5  planar  data.  However  Mack  noted  “the  major 
difficulty  in  the  use  of  the  forcing  theory  is  that  forced  disturbances  are  distinct  from  free 
disturbances,  and  the  process  by  which  the  former  becomes  the  latter  is  unknown”.  This 
coupling  between  freestream  (external)  sound  and  normal  modes  of  planar  supersonic 
boimdary  layer  has  been  theoretically  studied  in  [32],  [35].  It  was  found  that  the  first  and 
second  modes  are  synchronized  with  acoustic  waves  near  the  leading  edge  of  a  flat  plate,  i.e. 
frequency  and  wavenumber  of  the  normal  mode  coincide  with  those  of  an  external  acoustic 
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wave.  That  causes  a  strong  excitation  of  the  normal  modes  by  freestream  noise  through  the 
mechanism,  which  is  quahtatively  different  from  the  subsonic  case  studied  in  [13],  [18], 

(iii)  For  low  supersonic  Mach  numbers,  the  dominant  boundary-layer  disturbances  have 
relatively  low  frequency.  The  wind  tunnel  freestream  has  intensive  acoustic  disturbances  at 
the  corresponding  frequency  band  and  provides  large  initial  amphtudes  for  normal  modes  at 
the  instability  growth  onset  associated  with  the  lower  neutral  branch.  The  major  source  of 
these  disturbances  is  due  to  radiation  of  acoustic  waves  by  the  turbulent  boundary  layer  on 
the  wind  tunnel  nozzle  walls.  That  causes  transition  Reynolds  numbers  in  conventional 
supersonic  wind  tunnels  to  be  significantly  lower  than  ones  in  flight  conditions.  Hypersonic 
transition  wind  tunnel  testing  presents  a  different  situation  [26].  If  the  dominant  boundary- 
layer  disturbances  belong  to  the  second  mode,  they  have  relatively  high  frequency.  Since 
most  of  acoustic  energy  is  concentrated  in  the  low  frequency  band,  the  noise  intensity 
rapidly  decreases  with  increasing  frequency  until  the  signal  is  lost  in  the  instrumentation 
noise.  Therefore,  the  “noisy”  hypersonic  wind  tunnel  can  be  “quiet”  for  the  second  mode 
frequency  band.  In  this  situation,  freestream  turbulence  (vortical  disturbances)  and  entropy 
spottiness  (tenq)erature  disturbances)  may  play  inoportant  role  in  hypersonic  boundary-layer 
receptivity.  Mechanism  of  the  boundary-layer  mode  excitation  by  vorticity/entropy  waves  is 
discussed  in  Section  5. 

(iv)  Transition  process  is  described  by  the  initial  boxmdary  value  problem,  wfrich  requires 
distributions  of  the  normal  mode  anq)Htudes  at  the  instability  growth  onset.  These  initial  data 
are  determined  from  receptivity  analyses  conq)rising  the  identification  of  external 
disturbances,  the  determination  of  the  most  sensitive  regions,  and  prediction  of  the  normal 
mode  propagation  from  these  regions  to  the  instability  growth  onset.  All  these  steps  relate  to 
the  disturbance  field  in  the  region  located  upstream  from  the  lower  neutral  branch. 
However,  previous  studies  have  been  focused  on  stabihty  analysis  related  to  the  downstream 
region  only.  In  this  report,  we  show  that  the  disturbance  pre-history  of  hypersonic  bormdary 
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layer  is  more  complicated  than  in  the  moderate  supersonic  case  due  to  the  presence  of  the 
disturbance  spectrum  singularities. 

(v)  Recent  experimental  studies  [36]  relevant  to  the  non-linear  phase  of  supersonic  flat-plate 
transition  showed  that  boimdary-layer  disturbances  of  discrete  spectrum  are  able  to  generate 
intensive  acoustic  waves  of  continuous  spectrum.  Release  of  the  boundary-layer  disturbance 
energy  to  the  external  flow  can  affect  the  transition  onset  process.  This  mechanism  seems  to 
be  stronger  at  higher  Mach  numbers.  Since  the  PSE  analysis  does  not  accoimt  for  the 
continuous  spectrum  disturbances,  these  effects  cannot  be  predicted  using  the  PSE  code. 
Similar  comment  can  be  made  regarding  receptivity  problem  coupling  external  disturbance 
of  continuous  spectrum  with  the  boundary-layer  modes  of  discrete  spectrum.  To  capture 
these  mechanisms,  the  PSE  method  should  be  modified  or  combined  with  alternative 
methods. 

These  examples  show  that  the  extension  of  subsonic  and  low  supersonic  stability  concepts 
and  transition  prediction  methodology  to  hypersonic  speeds  is  not  a  trivial  problem 
Hypersonic  boundary  layer  receptivity,  instability  and  non-linear  mechanisms  reveal  new 
features  which  should  be  studied  independently.  A  multiple-mode  method  rather  than  a 
single-mode  approach  can  help  to  explain  mentioned  above  singular  featmes  of  hypersonic 
transition.  This  method  can  be  able  to  evaluate  receptivity  of  hypersonic  boundary  layers  and 
provide  initial  data  for  the  PSE  calculations.  That  will  help  to  realize  the  amplitude  method 
for  the  hypersonic  transition  prediction. 

This  report  addresses  the  initial  phase  of  the  hypersonic  boundary-layer  transition 
comprising  excitation  of  unstable  normal  modes  (receptivity  problem)  and  their  downstream 
evolution  fi-om  receptivity  regions  to  the  instability  growth  onset  (instabihty  pre-history 
problem).  In  Section  2,  we  formulate  basic  formahsm  of  the  multiple-mode  method 
following  to  previous  studies  relevant  to  subsonic  and  moderate  supersonic  speeds  [12], 
[37],  [38].  In  Section  3,  we  analyze  the  spectrum  topology  for  hypersonic  boundary  layer 
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including  the  continuous  spectrum.  Then  we  consider  normal  mode  characteristics  of  the 
hypersonic  boundary  layer  on  a  sharp  cone  and  identify  singular  regions  relevant  to 
synchronism  between  normal  modes  and  external  acoustic  waves,  synchronism  between  the 
first  mode  and  external  vorticity/entropy  waves,  and  synchronism  between  the  first  and 
second  modes.  In  Section  4,  we  study  coupling  between  the  first  and  second  mode  in 
correspondent  synchronism  region  and  estabhsh  the  inter-mode  exchange  rule.  In  Section  5, 
we  consider  synchronism  between  the  houndary-layer  modes  and  the  entropy/vorticity 
waves  and  estimate  receptivity  to  the  fireestream  turbulence  and  entropy  spottiness.  In 
Section  6,  we  summarize  results,  formulate  possible  scenarios  for  the  initial  phase  of 
hypersonic  transition  in  conical  and  planar  boundary  layers,  and  discuss  their  consistency 
with  experimental  evidences. 


2.  Basic  formalism 

In  this  section,  we  describe  the  basic  formalism  of  the  multiple-mode  method,  which 
provides  coupling  between  external  disturbances  and  boundary-layer  normal  modes 
(receptivity  problem)  as  well  as  coupling  between  difiFerent  normal  modes  (inter-mode 
exchange  problem).  In  the  fi’amework  of  linear  stability  theory,  receptivity  and  iater-mode 
exchange  mechanisms  are  due  to  non-parallel  effects  of  basic  flow.  The  multiple-mode 
analysis  is  based  on  the  following  formaUsm  developed  in  [12],  [37]-[39]. 

We  consider  a  two-dimensional  boundary  layer  mean  flow  with  the  velocity  components 
(U,V)  =  (U*,V*)IU*  ,  pressure  P=P*/(pluf)  and  temperature  T=T*IT*.  Due  to 
downstream  growth  of  the  boundary  layer  thickness,  the  mean  flow  profiles  are  fimctions  of 
the  longitudinal  coordinate  x  =  x*  I L* ,  i.e. 

U  =  U(x,y),  V  =  eV(x,y),  P  =  P(x),  T^T(x,y)  (2.1) 
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where  y  =  y*  !  5  is  coordinate  normal  to  the  wall  surface  ,  5  =  /  U*  is  bovmdary 

layer  scale,  and  Z*  is  typical  streamwise  distance  from  the  body  nose  to  the  flow  region  to 
be  analyzed.  In  many  cases,  this  distance  is  associated  with  the  lower  neutral  branch  locus. 


Unsteady  disturbance  field  is  represented  as  a  vector 


m/  t\  (  ^  dd  du  dv  dd  r 

^{x,y,t)  =  {u,  — ,  v,p,  e,  — ,  — ,  — ,  — ) 

ac  ^  (3c 


(2.2) 


where  t  =  t*Ul  1 5  \s  time;  u  and  v  are  streamwise  and  vertical  velocity  components;  p  is 
pressure  and  6  is  temperature.  For  harmonic  disturbances  of  a  given  angular  frequency 
CO  =  CO* 5 lU]  and  small  amplitude  r ,  the  disturbance  vector  is  expressed  in  the  form 

'F  =  r^(x,y)exp(-/fflt)  (2.3) 

The  disturbance  anq)litude  .4  is  a  solution  of  linearized  Navier-Stokes  equations,  which  can 
be  written  in  the  matrix-operator  form 

^  /M 

H^A  +  sH^  —  =  Q  (2.4) 

3c 

where  the  operator  is  expressed  as 

H,  =—L,  —  +  L—-H,-sH, 


Elements  of  the  9x9  matrices  4,  Zj,  and  are  fimctions  on  mean  flow 

characteristics.  Their  explicit  form  is  given  in  [12]  and  [38].  Elements  of  the  matrix  are 
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proportional  to  terms - ,  — ,  V  and  — .  Small  parameter  s  =  l/  characterizes  non- 

dx.  dx.  dy 

parallel  effects,  where  the  Reynolds  number  is  =  yju^L*  /  v*  .  If  the  disturbance 
amphtude  is  identified  in  a  certain  cross-section  x  =  x,,  and  disturbances  satisfy  no-shp 
conditions  on  the  wall,  then  the  boundary  conditions  are 


^(^0  7T)=  A(t) 

(2.5a) 

a4j(x,0)  =  243(x,0)  =  245(x,0)  =  0 

(2.5b) 

|i(x,>;)|<oo, 

(2.5c) 

Solution  of  Eq.  (2.4)  with  the  boundary  conditions  (2.5)  can  be  represented  as  an  expansion 
in  an  eigenvector  system  of  the  locally  homogeneous  problem 


r 

+  L  — 


=  H,A  +  eH^ 


(2.6a) 


4(x,0)  =  4(x,0)  =  ^5(x,0)  =  0 


(2.6b) 


|i(x,j)|<oo, 


(2.6c) 


The  biorthogonal  eigenvector  system  {/4^(x,7),  .6^(x,>;)}  is  defined  as  a  family  of 
solutions  for  the  following  problems  [38],  [39]: 


Direct  problem: 


\  Sy  j 


(2.7a) 

(2.7b) 

(2.7c) 


Adjoint  problem: 
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m.  ^ 


dy 


=  5„4(x,0)  =  5„6(x,0)  =  0 


(2.8a) 

(2.8b) 

(2.8c) 


Here  asterisk  denotes  a  conjugate  matrix  and  upper  bar  denotes  cottq)lex  conjugate  values. 
Eigenfunction  A^(x,y)  characterizes  an  amplitude  distribution  of  the  elementary  wave 


JT 

\ (>;.>’)exp[?’.<Jc> -lal],  if’Ax)  =  is'' J a(x)dx 


(2.9) 


The  eigenvectors  satisfy  the  following  orthogonahty  conditions 


(2.10) 


where  the  scalar  product  is  defined  as 


oo  9 

{A,B)=  ]im\e-‘^(A,B)dy,  =  ?c>  0  (2.11) 

Here  is  Kronecker  symbol,  if  one  of  the  eigenvalues  a,  belongs  to  the  discrete 
spectrum;  A^  =S{a-P)  is  delta-function,  if  both  eigenvalues  belong  to  the  continuous 
spectrum. 

Solution  of  the  problem  (2.4),  (2.5)  is  represented  as  the  eigenmode  decomposition 
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(2.12) 


a 


where  ^  denotes  summation  over  discrete  spectrum  and  integration  over  continuous 

a 

spectrum.  Substituting  (2.12)  into  Eq.  (2.4),  multiplying  by  Bp  and  accounting  for  the 

orthogonality  condition  (2. 10),  we  obtain  the  ordinary  diSerential  equation  system  for  the 
amplitude  coefficients  (x) 


=  (2.13) 

c„  =  at  x  =  Xo 

where  matrix  elements  are  determined  as 


(2.14) 


For  an  arbitrary  normalization  of  the  eigenfunctions  A^,B^,  these  matrix  elements  are 
e?q)ressed  as 


{hAA) 


(2.15) 


Diagonal  elements  relate  to  a  single  mode  evolution  in  the  non-parallel  boundary  layer. 
Non-parallel  stability  analyses  [12],  [41]  and  the  PSE  method  [19],  [20]  address  this  case. 
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Non-diagonal  elements  are  responsible  for  exchange  between  modes  with  eigenvalues 
a  and  .  If  one  of  these  modes  belongs  to  the  continuous  spectrum,  then  the  elements 

provide  coupling  between  external  disturbances  and  the  boundary-layer  modes  (receptivity 
problem).  If  both  modes  belong  to  the  discrete  ^ectrum,  these  elements  are  responsible  for 
the  exchange  between  normal  modes  (the  inter-mode  exchange  problem).  If  the  mean  flow  is 
parallel,  then  all  matrix  elements  =  0 ,  and  no  coupling  between  modes. 

In  the  framework  of  this  formalism,  the  receptivity  analysis  and  the  inter-mode  exchange 
analysis  comprises  the  following  steps: 

•  Identify  the  biorthogonal  eigenvector  system  {A^(x,y),  B^{x,y))  including  discrete  and 
continuous  spectrums  of  a(x,o)) . 

•  Specify  initial  conditions  at  x  =  . 

•  Simplify  the  equation  system  (2. 13)  neglecting  eigenmodes,  which  are  not  important  for 
particular  problem  (cutoff  procedure). 

•  Calculate  eigenvalues,  eigenfunctions  and  matrix  elements  (2. 14). 

•  Solve  the  cutoff  equation  system  for  the  anq)Utude  coefficients  (x) . 


3.  Disturbance  spectrum 


Analysis 

The  local  homogeneous  problem  (2.7)  can  be  reduced  to  the  well  known  linear  stabiUty 
problem  [12],  [38],  [39] 


dz„  _  ^ 

(3.1a) 

=-^a3  =Z„5  =0,  >;  =  0 

(3.1b) 

(3.1c) 
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where  ^a4->  ■^asf  ^  311^  G,,  is  6x6  matrix  depending  on  basic 

flow  characteristics,  fi'equency  co ,  eigenvalue  a ,  and  Reynolds  number  R  =  -yjulx*  /  v*  . 

Its  non-zero  elements  are  given  in  [41],  System  (3.1)  has  six  independent  solutions 
,  ^2  »•  •  •  5  ^6  "'^th  the  following  asymptotic  forms 


qj  at  y^co,  J  =  1,...,6 

where  Aj  are  solutions  of  the  characteristic  equation 

det||G(,  -  AE^  =  0  as  y  — >  oo 
Eq.  (3.3)  can  be  written  in  the  expUcit  form 

(6, .  - )[( - 1=  -  A" )  -  ]  =  0  , 

4„  =  Gl\  =  G^Cf  +Gt>G^  +G«Gr. 
*,.=G«Gf +G;>G”+G;‘Gf. 
*,2=Gr,»„=Gf 

Eq.  (3.4)  has  the  following  solutions 


(3.2) 


(3.3) 


(3.4) 


^11  ?  ^2,3 


(3.5) 


Speciiying  branches  as  Re(A,,  X^,  /I,)  <0,  we  can  express  the  asymptotic  forms  (3.2)  as 


?u  ^  .  ?,.4  ->  ,  ?«  ^  at  ^  00 


(3.6) 
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If  Re(A^.)5iO  for  all  y  =  l,  2, 3,  then  solution  of  the  problem  (3.1)  is  the  linear 
decomposition  of  three  vectors 


—Cyq^  +C^q^ 


(3.7) 


where  constants  Q ,  Cj  and  Cj  are  determined  from  the  wall  boundary  condition  (3.  lb) 


^1^11  ^3^13  ^5^15  ~  ®  (^•^) 

^1^31  '*’^3^33  "^^5^35  ~  ^ 

^1^51  "^^3^53  "^^5^55  ~  ^ 


Equation  system  (3.8)  has  non-trivial  solutions,  if  its  determinant  is  zero;  i.e. 


<ln 

^13 

^15 

^31 

^33 

^35 

=  0 

(3.9) 

^51 

^53 

^55 

The  dispersion  relationship  (3.9)  gives  a  discrete  spectrum  a„ico,x)  related  to  the  boundary 
layer  normal  modes. 

If  one  of  the  eigenvalues  Xj  has  zero  real  part,Re(Aj.)  =  0,  then  solution  of  the  problem 

(3.1)  is  a  linear  decomposition  of  four  vectors  and  describes  waves  of  continuous  spectrum. 
Analysis  [39]  shows  that  supersonic  boundary  layer  has  seven  branches  of  continuous 
spectrum:  three  branches  correspond  to  waves  propagating  upstream,  two  branches 
correspond  to  acoustic  waves  propagating  downstream  and  two  branches  correspond  to 
vorticity  and  entropy  waves  propagating  downstream  The  continuous  spectrum  branches 
and  discrete  eigenvalues  are  schematically  shown  in  Fig.  1.  Further  analysis  addresses  the 
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continuous  spectrum  waves  propagating  downstream  only.  Their  eigenvalues  belong  to  the 
upper  half  of  complex  a  -plane. 

If  =-k^  and  A:  >  0  is  real  parameter,  then  the  wavenumber  a^(co,k)  corresponds  to 
vorticity  waves  of  the  amplitude 

=  Ci^i  +C2^2  +^^3  (311) 

which  has  the  asyn:q)totic  form 

at  oo  (3. 12) 

For  large  Reynolds  numbers  R  =  ■^U*x*  /  vl  =  R^^yfx ,  the  vorticity  wavenumber  is 
approximated  as 


1  2  2 

a,=(o  +  i  -~  +  0{R-^)  (3.13) 

R 

If  /Ij  =  3nd  A:  >  0,  then  the  wavenumber  a^(cD,k)  corresponds  to  entropy  waves  of 
the  amplitude 


+^*3^3  +C^q^ 


(3.14) 


which  has  the  asyn^totic  form 


+€3^3^  +C,ql%^^  aty^oo  (3.15) 
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The  entropy  wavenumber  is  approximated  as 


a^=(o+i - - —  +  0{R^)  (3.16) 

(jR 

where  a  is  Prandtl  number.  According  to  Eqs.  (3.13)  and  (3.16)  vorticity  and  entropy 
waves  propagate  downstream  with  the  phase  speed  c^=o)l and  slowly 
dissipate  due  to  viscosity. 

Equation,  at  k>0,  gives  three  branches.  Two  of  them  correspond  to  acoustic 

waves  of  the  wavenumber  a^{(o,k)  propagating  downstream  in  a  supersonic  flow.  Their 
amplitude 


(3.17) 


has  the  asymptotic  form 

+  C,qf'e‘’  +  at  ^  »  (3. 18) 

For  the  Mach  munber  >  1  and  k  «R,  acoustic  waves  are  almost  neutral.  Their 
wavenumbers  belong  to  the  regions 

Re(o.,)Sr^  +  0(«");  Re(a„)>-^  +  0(ft-)  (3.19) 

As  k  increases,  their  decrement  Im(a„)  grows.  When  the  parameter  k  is  of  the  order  of 

the  Reynolds  number  R,  acoustic  waves  transform  into  pressure  waves,  which  rapidly 
dissipate  downstream  (see  Fig.  1). 
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Numerical  results 


Two-dimensional  boundary-layer  profiles  were  calculated  using  self- similar  solutions  of  the 
boundary-layer  equations  in  the  local  similarity  approximation  [42].  The  conical  mean  flow 
characteristics  were  obtained  using  the  Mangier  transformation  according  to  which  the  flat- 
plate  boundary-layer  thickness  dp  is  larger  than  the  conical  boimdary-layer  thickness  5^  by 

-\/3  .  In  all  calculations,  it  is  assumed  that  the  fluid  is  a  perfect  gas  with  a  constant  specific 
heats  ratio  y  =  1-4  and  Ihandtl  number  cr  =  0.72 .  The  viscosity  temperature  dependency  is 
approximated  by  Sutherland’s  law 


(1  +  ^) 
(r+5) 


where  iS'  =  110/r/  for  air  temperature  measured  in  degrees  Kelvin;  //  =  //*///*  is 
nondimensional  viscosity.  Mean  flow  calculations  were  performed  for  the  conditions  related 
to  the  AEDC  Tunnel  B  stability  experiments  on  a  7  degree  half  angle  sharp  cone  at  zero 
angle  of  attack  [25].  Comparison  between  theoretical  and  experimental  distributions  5 ^{x) 
are  shown  in  Fig.  2  for  the  local  Mach  number  =  6.8 ,  local  ten^eratine  =  70.26  K; 
the  cone  surface  teirqrerature  is  (the  insulated  wall  condition).  It  is  seen  that  the 

local  self-similar  solution  well  correlates  experimental  data.  Consequently,  the  viscous- 
inviscid  interaction  and  the  non-uniform  distribution  of  the  surface  temperature  near  the 
cone  tip  can  he  neglected  in  this  case. 


Integration  of  the  stability  equations  (3.1)  is  performed  using  the  numerical  procedure 
similar  to  that  developed  by  Mack  [21].  The  analytic  external  solutions  (3.7)  provide  the 
initial  conditions  for  the  linearly-independent  solutions  at  the  upper  boimdary  T]=rj^,  where 

Ti  =  y  /  ^fx  is  self- similar  boundary-layer  variable.  The  system  (3.1)  is  integrated  fi^om  the 
upper  boundary  77^  to  the  wall  surface  tj  =  0  using  the  fourth-order  Runge-Rutta  method 
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and  the  orthogonalization  procedure.  The  eigenvalue  search  is  performed  by  the  Newton 
iteration  method.  Calculations  are  conducted  for  the  planar  boundary  layer.  These  results 

have  been  adjusted  by  Vs  for  the  comparison  with  the  conical  boundary-layer  stability  data 
[25].  That  is,  the  eigenvalues  a  and  the  Reynolds  number  R  have  been  increased  by  Vs  . 

Eigenmode  characteristics  and  boimdaries  of  continues  spectrum  are  shown  in  Figures  3a 
and  3b  for  the  local  Mach  number  -  6.8,  local  temperatme  =  1026  K,  the  cone 

surface  temperature  T^  =  T^,  and  frequency  parameter  =  10^ .  These  parameters 
correspond  to  stability  measurements  [25]  performed  in  the  AEDC  Tunnel  B  on  a  7  degree 
half  angle  sharp  cone  at  zero  angle  of  attack.  The  &st  and  second  mode  eigenvalues  tend  to 
the  lower  and  upper  boundaries  of  acoustic  spectrum  respectively,  when  the  Reynolds 
number  is  R<  800  (Region  1  in  Fig.  3  a).  Asymptotic  analysis  [32]  and  numerical 
calculations  conducted  in  a  wide  range  of  Mach  niuubers  and  wall  temperatures  showed  that 
this  upstream  asyn^totic  behavior  is  typical  for  hypersonic  boundary  layers  and  can  be  used 
for  classification  of  normal  modes.  We  will  call  the  normal  mode  with  the  phase  speed 
— >  1  + 1  /  at  X  — >  0  as  the  first  mode,  and  the  normal  mode  with  the  phase  speed 

->  1  - 1  /  at  X  —>  0  as  the  second  mode.  Note  that  this  classification  is  different  from 

the  Mack’s  classification.  To  avoid  misunderstanding  we  will  call  the  first/second  mode  of 
Mack’s  classification  as  the  Mack’s  first/second  mode. 

Analysis  [32]  showed  that  synchronization  of  the  first  and  second  mode  waves  with  acoustic 
waves  of  the  phase  speed  =  1  ±  1  /  can  cause  a  strong  excitation  of  instability  near  the 
nose  region.  In  conventional  wind  tunnels,  sound  is  generated  by  vortices  propagating 
downstream  in  the  turbulent  boundary  layer  on  the  nozzle  and  test  section  walls.  Since 
acoustic  waves  from  these  moving  sources  have  the  phase  speed  c  <  1  - 1  /  ,  they  do  not 

satisfy  the  synchronism  conditions.  However,  their  interaction  with  the  nose  tip  induces  a 
local  source,  which  generates  “dangerous”  disturbances  of  the  phase  speed  c  =  1  ±  1  /  ,  as 

schematically  shown  in  Fig.  4.  For  a  planar  configuration  (flat  plate  or  hollow  cylinder),  this 
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local  source  has  cylindrical  shape.  For  the  conical  flow,  it  is  concentrated  near  the  cone  tip. 
Since  the  point  source  on  the  cone  tip  is  weaker  than  the  line  somce  on  a  plate  leading  edge 
and  downstream  attenuation  of  acoustic  waves  from  the  point  source  is  stronger  than  that 
from  a  hne  source,  receptivity  of  a  planar  configuration  to  acoustic  disturbances  should  be 
higher  than  that  of  a  conical  configuration.  These  arguments  are  consistent  with 
experimental  observations  [33]. 

Near  the  point  « 1515  (Region  2  in  Fig.  3a)  the  first  mode  is  synchronized  with  vorticity 
and  entropy  waves  of  the  phase  speed  =  1.  In  this  region,  the  freestream  vorticity  and 

entropy  disturbances  may  effectively  generate  the  first  mode  and  conversely  the  first  mode 
may  induce  the  vorticity  and  entropy  disturbances.  These  mechanisms  are  discussed  in 
Section  5.  The  first  mode  mass  flux  distribution  Qi(y,R)  is  shown  in  Fig.  5.  Near  the 
synchronism  point,  jR«1515,  oscillating  “tail”  indicates  that  the  first  mode  penetrates 
outside  the  boundary  layer.  This  feature  can  be  exploited  to  identify  the  synchronism  region 
2  by  measiuing  correlation  between  hot-wire  signals  inside  and  outside  the  botmdary  layer. 

In  the  vicinity  of  the  point  R  «  1750  (Region  3  in  Fig.  3a),  phase  speeds  of  the  first  and 
second  modes  are  very  close  to  each  other.  As  will  be  shown  in  Section  4,  the  discrete 
spectrum  has  two  branch  points  in  this  region.  Near  the  branch  points  a  strong  exchange 
between  the  boundaiy-layer  modes  can  be  caused  by  the  mean  flow  nonparallel  effect. 
Mechanism  of  this  exchange  will  be  discussed  in  Section  4.  The  synchronism  is  associated 
with  local  maximum  of  the  first-mode  mass  flux  distribution  shown  in  Fig.  5. 

The  spectrum  topology  may  be  different  from  that  shown  in  Figures  3a  and  3b.  As  an 
example.  Figures  6a  and  6b  illustrate  the  eigenvalue  behavior  for  the  local  Mach  number 
=  5.5 ,  wall  teroperature  ratio  7^  /  7^  =  0.1  and  disturbance  frequency  F  =  10"^ . 
Downstream  the  branch  point,  «  3600  (Region  3),  the  first  mode  gets  unstable  whereas 
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the  second  mode  gets  more  stable.  In  Section  4,  we  study  the  spectrum  behavior  near  the 
branch  points  and  discuss  cases  corresponding  to  jumps  from  one  topology  to  another. 


4.  Coupling  between  the  first  and  second  modes 

Analysis 

Figures  3a  and  6a  show  that  phase  speeds  of  the  first  and  second  modes  are  very  close  to 
each  other  in  the  synchronism  region  3.  The  distmbance  spectrum  topology  in  this  region 
has  been  studied  by  Guschin  and  Fedorov  [30],  [3 1].  It  was  found  that  a  quaUtative  local 
behavior  of  eigenvalues  is  described  by  the  following  single  dispersion  relation 

=  (4.1) 

where  a  =  a- and  a  =co-co^  =  F(R -R^)  is  local  eigenvalue  and  angular  frequency 
respectively;  and  are  constants.  The  eigenvalues  of  ^  given  frequency  F 

are  schematically  shown  in  Fig.  7  for  the  case  of  <0,  a^>a^ >0.  Solution  of  Eq.  (4. 1)  has 
two  branches 


^1,2  ~ 


a^  +^2 


FR+F(a,-a,fF^R^+a, 


The  branch  points  are  defined  as 


Fia^-a^) 


(4.2) 


(4.3) 


They  are  located  on  the  real  axis  of  complex  R  -plane.  If  the  branch  points  are  bypassed 
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from  the  upper  side  as  shown  in  Fig.  7,  then  the  first  mode  is  stable,  Im(aj )  >  0 ,  and  the 
second  mode  is  unstable,  Im(a2 )  <  0,  in  the  region  <  i?  <  .  This  spectrum  topology 

is  consistent  with  the  numerical  result  illustrated  by  Figures  3  a  and  3b.  As  will  be  shown 
below,  the  branch  point  is  slightly  shifted  below  the  real  axis  of  the  complex  -plane 
due  to  viscous  effects.  As  the  Reynolds  number  increases  along  the  real  axis,  this  branch 
point  is  bypassed  from  the  upper  side  similar  to  the  case  shown  in  Fig.  7.  Figures  6a  and  6b 
illustrate  another  case  when  the  branch  point  is  shifted  to  the  upper  half  of  the  complex 
-plane  and  bypassed  from  the  lower  side.  The  first-mode  wave  gets  unstable  downstream 
from  the  branch  point  ,  whereas  the  second  mode  gets  stable. 

Now  we  consider  disturbance  behavior  near  the  branch  point  Xjj  =  !  Rj^Y .  Eigenvalues 

of  the  first  and  second  modes  can  be  approximated  as 

=  ^0  ±  aVx+...  (4.4) 

where  A"  =  x  -  Xjj  is  local  variable.  Using  the  basic  formahsm  of  Section  3,  we  can  describe 
the  coupling  between  the  first  and  second  mode  in  the  two-mode  approximation.  I.e. 
solution  of  Eq.  (2.4)  is  approximated  as 

i  =  Cj(x)ij(x,y)exp[^j(x)]-i-C2(x)i2(^,J^)exp[^2(x)]  (4.5) 

X 

^1,2  (4.6) 

^b\ 


where  the  eigenfrmctions  solutions  of  the  homogeneous  problem 


=0 


(4.7a) 
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Ajj(x,0)  =  43(x,0)  =  45(x,0)  =  0 
1 4  0,y^oo;j  =  l,2 


(4.7b) 

(4.7c) 


The  ampUtude  coefficients  Cj  2  are  solutions  of  the  ordinary  differential  equation  system 


(4.8a) 
(4.8b) 

where  the  matrix  elements  are  defined  as 


aCj 

cbc 

dc^ 

dx 


=  cJV,,  +C2lTi2exp(^2  -<Pi) 


=  €2^22  +c,ir2iexp(^i  -^2) 


W  = 


{hJj.B,} 


,  J  —  I,  2,,  k  —  1, 2 


(4.8c) 


In  the  vicinity  of  the  branch  point,  Jif  =  0,  the  operators  and  eigenfimctions 

are  expanded  in  the  following  series 


•^0  -  ^00  +-^oi+"-  (4.9a) 

H2  =  + XH21+...  (4.9b) 

Aj  =  +(«j.  -  ao)A^^^+...  ;  J  =  1,2  (4.9c) 

Substituting  (4.4)  and  (4.9)  into  Eqs.  (4.7)  we  obtain  the  following  problems 

(H,,+ia,H2o)A^°^  =0  (4.10a) 

4”^  (:>c,0)  =  4“^  (^,0)  =  ^  (^,0)  =  0  (4. 10b) 
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1  .y  ^ 

(4.10c) 

(4.11a) 

^«(x,0)  =  4'^(x,0)  =  ^f>(x,0)  =  0 

(4.11b) 

|i<‘^(x,j);)|^0,  j^oo 

(4.11c) 

(4.12a) 

A^^\x,0)  =  ^3^^(x,0)  =  Aj^\x,0)  =  0 

(4.12b) 

lA^^^(x,y)j^0,y^oo 

(4.12c) 

A  non-trivial  solution  of  the  problems  (4.11)  and  (4.12)  exists,  when  the  right-hand  side  of 
Eqs.  (4.11a)  and  (4.12a)  is  orthogonal  to  the  adjoint  problem  solution  5^”^ ;  i.e. 


(4.13) 

(4.14) 


Eq.  (4.13)  can  be  used  for  the  branch  point  search.  In  turn,  the  constant  A  can  be 
determined  from  Eq.  (4. 14). 


Similar  e?q)ansions  can  he  made  for  the  adjoint  problem  solution  B .  Substituting  them  into 
Eq.  (4.8c)  we  obtain  the  following  form  of  the  matrix  elements 


/_iy+A-i 

Wj,  =  ^  -  ;  j  =  1,2;  k  =  l,2 


(4.15) 
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In  the  inner  region  with  the  variable  =  s  =  0(1) ,  the  equation  system  (4.8)  is  written 
as 


C/Cj 


4^ 


-Cj  +  Cj  exp 


J_ 

4^ 


c,  exp| 


(4.16a) 

(4.16b) 


Since  this  system  does  not  depend  on  the  mean  flow  profiles,  its  solution  can  be  used  to 
model  the  inter-mode  exchange  for  broad  class  of  shear  layers  with  the  eigenvalue  behavior 
similar  to  that  given  by  Eq.  (4.4).  Equations  (4.16)  can  be  reduced  to  the  Airy  equation 
using  the  following  transformation 

/(4)  =  Ci(^)exp(A(^))  +  C2(^)exp(-A(^))  (4.17a) 

A(^)  =  |/A^''"  (4.17b) 


Substituting  (4. 17)  into  (4. 16)  we  obtain 

^  +  A^^  =  0  (4.18) 

Sunilar  equation  has  been  derived  by  Fedorov  and  Khokhlov  [32]  using  an  alternative 
asymptotic  analysis  of  Eq.  (2.4). 

Asymptotic  behavior  of  the  Airy  equation  solution  at  |^->  oo  depends  on  the  orientation  of 
Stokes  lines  /j ,  j  =  l,  2, 3,  which  are  defined  by  the  equation  Re[A(^)]  =  0  in  the  conq)lex 
^ -plane.  We  consider  the  Stokes  lines  topology  schematically  shown  in  Fig.  8.  This 
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topology  corresponds  to  the  case,  when  the  first  and  second  modes  are  almost  neutral  in  the 
upstream  region  <0.  Downstream  jfrom  the  branch  point  one  of  these  modes  gets 
essentially  unstable  and  another  gets  stable.  Numerical  exartq)les  of  such  spectrum  sphtting 
are  shown  in  Figures  3a,  3b  and  6a,  6b. 

Following  to  the  analysis  of  Fedoruk  [43]  we  introduce  the  canonical  domains 

D/Jj  e^j 

and  specify  branches  of  the  fimction  A(^)  so  that  Im(A)  >0  at  ^  elj  ^Dj.  Then  the  real 
part  Re(A)  is  positive  to  the  right  of  the  Stokes  line  Ij ,  and  it  is  negative  to  the  left  of  the 
Stokes  line  Ij .  Solution  of  the  Airy  equation  (4. 18)  has  the  asynq)totic  form 

/  [Cj  exp(A)  +  bj  exp(- A)] ,  1 41  ^  oo  (4.19) 

where  ^  eDj',  constants  Sj  satisfy  the  following  conditions 


\Sj\=  1,  and  arg(A^.^  =  0  st  ^  el j-,  j  =  1,2,3 
Constants  Uj  and  bj  are  coupled  by  the  relations  [43] 

(4.20a) 
(4.20b) 


at  j  =  \,2 


’o  r 

a. 

1 

A. 

=  exp(-/;r  /  6) 

1  / 

3 

A. 

j+i , 


=  exp(-/;r  /  6) 


0  1 
1  i 
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Eqs.  (4.20)  deteramie  the  exchange  between  the  first  and  second  modes  as  the  longitudinal 
coordinate  x  passes  through  the  interaction  region  |x-Xjj|=  and  intersects  one  of 

the  Stokes  lines.  Since  the  transition  matrix  has  a  universal  form,  Eq.  (4.20)  establishes  the 
inter-mode  exchange  rule  for  various  weakly  nonparaUel  mean  flows,  which  have  similar 
topology  of  the  discrete  spectrum  This  rule  can  be  expressed  in  more  generic  form  as 


^stable 

umtcble 


exp(-/^) 

o 


0  1 
1  i 


stable 


^unstable 


(4.21) 


where  (  )  and  ( A*^,^ ,  )  are  amplitudes  of  stable  and  unstable  mode 

upstream  and  downstream  fi’om  the  interaction  region  respectively.  For  the  eigenvalues 
shown  in  Figures  3a  and  3b,  the  ampHtudes  A~^,^ ,  A^^,^  correspond  to  the  first  mode,  and 
the  ampHtudes  correspond  to  the  second  mode.  For  the  case  shown  in 

Figures  6a  and  6b,  the  ampHtudes  A~^,^ ,  A^^,^  correspond  to  the  second  mode,  and  the 
ampHtudes  correspond  to  the  first  mode.  The  inter-mode  exchange  diagram 

shown  in  Fig.  9  indicates  that  the  stable-mode  wave  excites  the  imstable-mode  wave  of  the 
same  anq)Htude,  In  turn,  the  unstable  mode  excites  the  stable  mode  and 

propagates  downstream  with  the  same  initial  ampHtude.  In  both  cases,  the  unstable  mode 
occurs  in  the  downstream  region. 

In  the  first  order  approximation,  the  inter-mode  exchange  does  not  depend  on  the  small 
parameter  e ,  which  characterizes  nonparaUel  effects.  As  e  tends  to  zero,  the  interaction 
region  is  expanded  as  ,  whereas  the  local  disturbance  aitq)Htude  grows  as  s~^'^ .  Both 
trends  compensate  each  other  so  that  coupling  between  the  input  and  output  anq)Utudes 
does  not  change. 
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Numerical  results 


Numerical  analysis  of  the  dispersion  relation  (3.9)  is  extended  to  the  complex  values  of  the 
Reynolds  number  R .  For  self-similar  approximation  of  the  mean  boundary-layer  flow,  this 
extension  is  trivial,  because  the  mean  flow  profiles  do  not  depend  on  the  Reynolds  number 
and  the  stability  matrix  from  Eq.  (3.1a)  depends  on  the  Reynolds  number  in  parametrical 
explicit  form  only.  In  this  case  the  eigenvalue  problem  (3.1)  can  be  solved  for  complex 
values  of  R  using  the  same  numerical  procedures  as  in  cases  of  real  R .  The  branch-point 
Reynolds  number  is  calculated  Jfrom  Eq.  (4.13),  which  can  be  written  as 


(/Tj  =  ^  at  y  =  1, 2  and  (4.22) 

Eq.  (4.22)  is  solved  using  the  Newton  iteration  technique.  Distributions  Re[i^j(F)]  and 
Im[i^j(F’)]  are  shown  in  Figures  10a  and  10b  for  the  case:  local  Mach  number  =  6.8, 
wall  tenqrerature  ratio  =1.  It  is  seen  that  the  branch  point  is  located  in  the  lower 

half  of  the  complex  7? -plane  near  the  real  axis.  Since  the  branch  point  is  bypassed  fi"om  the 
upper  side,  the  spectrum  q)litting  is  similar  to  the  case  shown  in  Fig.  3a  and  3b. 

Figures  11  and  12  show  that  the  real  part  of  the  branch  point,  Re(ii^j) ,  is  close  to  the  low 
branch  of  the  second-mode  neutral  curve  predicted  by  Simen  and  Dallmann  [28]. 
Discrepancy  between  ejq)erimental  data  and  stability  calculations  [28]  observed  near  the 
lower  neutral  branch  may  be  due  to  singular  behavior  of  the  disturbance  field  in  the 
interaction  region,  R- R^^  =  OiR^^),  corresponding  to  the  slow  variable  range 
|3c-Xji|=  According  to  the  asymptotic  formula  (4.19),  the  disturbance  amplitude 

has  a  local  maximum  proportional  to  £~^'^  (see  also  the  eigenfunction  in  Fig.  5).  This  local 
increase  of  the  disturbance  an^Utude  can  be  misinterpreted  as  the  downstream  growth  of  the 
boundary-layer  instabihty.  Relatively  small  increments  calculated  fi’om  experimental 
measurements  near  the  lower  neutral  branch  may  relate  to  this  effect. 
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Figures  13a  and  13b  show  the  branch-point  Reynolds  number  as  a  function  of  local 
Mach  number  at  various  wall  tenq)erature  ratios  for  the  disturbance  frequency 

parameter  F  =  10^.  For  the  cooled  wall  case,  =  01,  the  imaginary  part  Im(/?jj) 

changes  sign  at  M^,  «  5.9 .  In  the  region,  ,  the  branch  point  belongs  to  the  lower 

half  of  the  conq)lex  R  -plane.  As  the  Reynolds  number  increases  along  the  real  axis,  the 
branch  point  is  bypassed  from  the  upper  side  as  schematically  shown  in  Fig.  7.  I.e.  the 

second  mode  is  stable  and  the  first  mode  is  xmstable  in  the  downstream  region  R  >  Re(i?jj ) . 
This  topology  corresponds  to  the  case  shown  in  Fig.  14a  and  14b  for  =  5.95  and 
!  '^ad-  ■  I®  tib®  region  ,  the  branch  point  belongs  to  the  upper  half  of  the 

corrq)lex  R  -plane,  and  it  is  bypassed  from  lower  side.  The  second  mode  is  unstable  and  the 
first  mode  is  stable  in  the  downstream  region  R  >  Re(/?j,j) .  This  topology  corresponds  to 
the  case  shown  in  Figures  6a  and  6b  for  =  5.5  and  !  T^=  0.1 .  Comparing  these  cases 
we  conclude  that  small  change  of  the  Mach  number  near  the  intersection  point  leads  to 

the  junq)  from  one  spectrum  topology  to  another.  However,  the  inter-mode  exchange  rule 
(4.21)  shows  that  the  imstable  mode  amphtude  is  approximately  the  same  in  both  cases;  i.e. 
the  initial  amphtude  of  instabihty  is  weakly  sensitive  to  this  jump 

If  the  real  axis  passes  through  the  local  region  |x-Xjj|=  i.e.  |Im(Xjj)|= 

and  |Im(/?j,j)|=  then  the  normal  mode  decon:q)osition  has  singular  behavior  in  the 

vicinity  of  the  branch  point  R^^.  According  to  Eq.  (4.19)  the  normal  wave  anq)htude  tends 
to  infinity  as  R-^  R^^.  Moreover,  the  group  velocity  of  normal  waves 

r^=Re 

rapidly  decreases  to  minus  infinity  and  then  jun^s  to  plus  infinity  as  the  Reynolds  number 
crosses  the  branch  point  from  upstream  to  downstream  region.  These  singularities  can  cause 


dco 

da 


1 


1 


Re(^a  /  dcD)  Fda,.  /  dR 


(4.23) 
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significant  difficulties  in  the  PSE  integration  near  the  branch  points  of  discrete  spectrum. 
Standard  stability  calculations  using  the  or  COSAL  code  can  also  be  affected,  because 
these  codes  determine  the  instability  trajectory  from  the  group  velocity  analysis. 

5.  Receptivity  to  vorticity/entropy  waves 

Problem  formulation 

In  Region  2  indicated  in  Fig.  3a  (see  also  Fig.  6a),  the  first-mode  eigenvalue  intersects 
branches  of  vorticity/entropy  waves  with  the  phase  speed  —  1.  The  intersection  region  is 
shown  on  an  enlarged  scale  in  Figures  15a  and  15b  for  the  local  Mach  number  =  6.8, 
wall  ten:q)erature  ratio  !  T^  =  l  and  disturbance  frequency  F  =  10^ .  Figure  16  shows 
the  spectrum  topology  in  the  complex  a  -plane.  Since  the  entropy-wave  branch  coincides 
with  the  vorticity-wave  branch  as  the  Reynolds  number  both  branches  can  be 

represented  in  the  complex  a  -plane  as  one  branch  cut.  The  first-mode  eigenvalue  a,  enters 
into  the  left  side  of  the  vorticity/entropy  branch  cut  at  the  Reynolds  number  -  0  and 

departs  from  its  right  side  at  the  Reynolds  number  R  =  R,+0,  where  «1515  is 

synchronism  point.  The  eigenvalue  imaginary  part  has  the  jurtq)  A(a,. )  «  0.34  •  10"^ ,  whereas 
its  real  part  is  smooth,  i.e.  =  ay^{R,  +0) .  Eigenfimctions  of  the  x  -component 

velocity,  ^(t/)  =  Re[z,(7)],  are  shown  in  Fig.  17  for  various  Reynolds  numbers  in  the 
vicinity  of  synchronism  point  R,.  These  fimctions  are  normalized  by  the  condition 
max|M(;7)|=  1.  It  is  seen  that  velocity  oscillations  penetrate  into  external  flow  as  the 
Reynolds  number  tends  to  the  synchronism  point  R,  (see  also  Fig.  5).  The  eigenfimction  at 
R  =  R-Q  (Line  3  in  Fig.  18)  is  essentially  different  from  the  eigenfimction  at  R=^R+Q 
(Line  4  in  Fig.  18).  This  is  due  to  the  fact,  that  asymptotic  behavior  of  the  first  vector  from 
(3.6),  at  y  -»c»,  is  very  sensitive  to  the  eigenvalue  junp  because 

A/lj  « -JiRAa  =0(1).  Consequently,  the  first  mode  should  be  treated  as  two  different 
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modes:  one  mode  corresponds  to  the  region  R<R^  and  another  to  the  region  R>  R,. 

VJq  assume  that  nonparallel  effects  can  cause  significant  interaction  between  the  first  mode 
and  vorticity/entropy  waves  in  the  synchronism  region  2  (see  Figures  3a  and  6a).  Neglecting 
other  modes  of  discrete  and  continuous  spectrum  we  represent  solution  (2. 12)  in  the  form 

oo 

A  =  c^(x)A,(x,y)Qxp((p^)  +  jc^(x,k)A^(x,y;k)exp((pJdk  + 

0 

oo 

+jcXx,k)R(x,y;k)Qxp(^^)dk  (5.1a) 

0 

XX  X 

<p^(x)  =  is~^  jai(x)dx,  q>^{x)  =  is~^  ja^(x,k)dx,  <p^(x)  =  is~^  j a^(x,k)dx  (5.1b) 

■*0  -*0  ^0 


where  subscripts  “v”  and  “e”  denote  vorticity  and  entropy  waves  respectively.  The 
equation  system  (2. 13)  is  written  as 

dc  f 

=  CjlTij  +jc^(x,k)W,^(x,k)exp[q)^(x,k)-(p^(x)yik  + 

CDC  ^ 

oo 

+\c^(x,k)Wi^(x,k)Qxp[^^{x,k)-^^(x)}Jk  (5.2a) 

0 

dc(k)  ■? 

exy[(pXx,k)  -  (P\(x)]  +  \cAx,k')W^(x,k,k')&ci)[(pAx,k')  -  (pAx,k)']fik'  + 
ax  ^ 

oo 

+\cXx,k')W^Xx,1i,k')Qxp[(p^{x,k')-(pXx,k)yik'  (5.2b) 

0 


dcSk) 

dx 


=  <^1^1  ^xy{(Peix,k)  -  ^7i(af)]  + J cAx,k’)W^Xx,k,k')Qxip[(pXx,k')  -  (pXx,k)W  + 

0 
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(5.2c) 


0 

Accounting  for  the  spectrum  topology  shown  in  Figures  15  and  16,  we  can  formulate  the 
following  \mit  problems  for  the  equation  system  (5.2): 

(i)  In  the  cross-section,  x  =  Xo,  located  upstream  from  the  synchronism  point,  x  =  x,, 
disturbance  spectrum  contains  only  the  first-mode  wave  of  a  unit  an^htude;  i.e. 

Cj(Xo)  =  l,  c„(x„,A:)  =  0  and  c^(Xo,A:)  =  0  forany  A:  e[0,oo)  (5.3) 

The  first-mode  wave  comes  into  the  synchronism  region  and  excites  vorticity/entropy  waves 
propagating  downstream.  Find  the  ampUtude  distributions  c^{x,k)  and  c^(x,k),  as  well  as 
the  fibrst-mode  wave  anq)htude  Cj(x)  in  the  downstream  region  x>x,. 

(ii) In  a  cross-section,  x  =  located  upstream  the  synchronism  point,  x  =  x, ,  disturbance 
spectrum  contains  only  one  vorticity  wave  of  the  wavenumber  =  a where  k, 
satisfies  the  synchronism  condition  (x, ,  A, )  =  a j  (x,  -I-  0) ;  i.e. 

=  c„(Xo,A)  =  ^(A-A.),  and  c^(Xj,A)  =  0  forany  A  €[0,c»)  (5.4) 

The  vorticity  wave  comes  into  the  synchronism  region  and  excites  the  first-mode  wave.  Find 
the  first-mode  wave  an^Utude  in  the  downstream  region  x  >  x, . 

(iii)  In  a  cross-section,  x  =  x„,  located  upstream  from  the  synchronism  point,  x  =  x,, 
disturbance  spectrum  contains  only  one  entropy  wave  of  the  wavenumber  =a^(A.), 
where  A.  satisfies  the  synchronism  condition  a^(x.,A.)  =  a^(x,  +0) ;  i.e. 
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Cj(Xo)  =  0,  cXy^^,k)  =  Sik-k,),  and  c^(Xo,A:)  =  0  forany  A:  e[0,oo) 


(5.5) 


The  entropy  wave  comes  into  the  synchronism  region  and  excites  the  JBrst-mode  wave.  Find 
the  first-mode  wave  amplitude  in  the  downstream  region  x  >  x, . 

Rough  estimates 

Due  to  the  presence  of  continuous  spectrum  it  is  difficult  to  solve  Eqs.  (5.2)  even  for 
relatively  simple  boundary  conditions  (5.3)-(5.5)  of  the  unit  problems  (i),  (ii)  and  (iii).  We 
guess  these  problems  can  be  simplified  using  asymptotic  analysis  of  the  matrix  elements 
W^{x,k,k’) ,  W^^(x,k,k') ,  W^^(x,k,k')  and  W^^(x,k,k')  as  well  as  integrals  over  the 

continuous  spectrum  branches  in  the  vicinity  of  synchronism  point.  This  could  be  a  subject 
for  future  research.  However,  rough  estimates  of  coupling  between  the  first  mode  and 
vorticity/entropy  waves  can  be  made  in  the  following  way.  We  introduce  a  virtual  upper 
boimdary,  y  =  y^,  located  in  the  inviscid  flow  region  between  the  upper  boundary-layer 
edge  and  the  shock.  On  this  boundary,  we  impose  the  homogeneous  boimdary  conditions 


2.1  =2„3  =  0,  y  =  y,  (5.6) 

Now  the  problem  (3.  la),  (3.  lb)  and  (5.6)  has  the  non-trivial  solution 

2.  ~  ^1^1  +^2^2  ■*■^4^4  ■*■^3^5 

if  the  determinant  specified  by  boundary  conditions  is  zero,  i.e. 
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The  dispersion  relation  (5.8)  leads  to  splitting  of  the  continuous  spectrum  branches  (see  Fig. 
1)  into  the  discrete  subsets.  As  the  upper  boxmdary  tends  to  infinity,  >>^->^00,  these  subsets 
thicken  back  into  correspondent  branches  of  continuous  spectrum. 

Numerical  analysis  of  the  problem  (3.  la),  (3.  lb)  and  (5.6)  shows  that  the  discrete  mode  with 
the  eigenvalue  a^,  which  belongs  to  the  vorticity/entropy  wave  subset  and  satisfy  the 
synchronism  condition,  cc^(R,  +0)  as  y,  ^  00,  transforms  into  the  first  mode  in 

the  region  R>  R,  .  In  a  similar  way,  the  first  mode  with  the  eigenvalue  coming  fi’om  the 
upstream  region,  R<R,,  transforms  into  the  mode  with  eigenvalue  which  belongs  to 
the  vorticity/entropy  wave  subset  and  satisfies  the  synchronism  condition 
a„(R,)->  a^(R,  -  0)  as  -4^  00 .  These  transformations  are  shown  in  Figures  18a  and  18b 

for  the  upper  boundary  locus  y  *  =  AS* ,  where  S*  is  boundary  layer  displacement  thickness. 
The  region  between  the  upper  boundary  layer  edge  y  =  y^  ,  and  the  virtual  boundary 
y  =  y^,  contaitts  about  12  wavelengths  of  the  vorticity/entropy  oscillations  in  the  y- 
direction,  i.e. 


(5.9) 

K 

where  the  parameter  k,  correspond  to  the  synchronism  condition.  Consequently,  this  case  is 
a  good  approximation  of  the  limit  y^  ->oo.  Li  the  upstream  region  R<R,  « 1515,  the 
imaginary  part  of  the  normal  mode  eigenvalue  (solid  line  in  Fig.  18b)  is  very  close  to  the 

imaginary  part  of  the  vorticity/entropy  wave  eigenvalues  (dashed  line),  which  satisfy  the 
synchronism  conditions 

av(^>*v.)  =  «i(^+9).  ««(^,^e.)  =  ai('^+0)  (5.10) 
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These  findings  lead  to  the  followmg  conclusions: 


•  In  the  rough  (of  the  order  of  magnitude)  approximation,  the  first  mode  coming  fi’om  the 

upstream  region,  transforms  into  vorticity  and  entropy  waves  with  the 

wavenumber  satisfying  the  synchronism  condition 

a^(K,K.)  =  a,iK-0),  a^(K,K,)  =  a,(Il-0)  (5.11) 

These  vorticity/entropy  waves  propagate  downstream  and  slowly  dissipate  due  to 
viscosity.  Since  this  mechanism  is  not  directly  coupled  with  excitation  of  unstable  normal 
modes,  it  plays  a  minor  role  in  transition  process. 

•  In  the  rough  approximation,  the  vorticity/entropy  waves  coming  fi’om  the  upstream 
region,  R<R,,  and  satisfying  the  synchronism  condition  (5.10),  transform  into  the  first 

mode  propagating  downstream,  R>  R,.  If  the  vorticity/entropy  eigenfimctions  A^{x,y) 
and  A^(x,y)  are  normalized  similar  to  the  first  mode  eigenfimction  A^(x,y),  say  the 
disturbance  pressure  an^litude  on  the  wall  A^^(x,0)  =  ^^^(^,0)  =  ^j4(x,0)  =  1,  then  the 
vorticity  wave  of  the  an:q)litude  c^(R,-0)  generates  the  first  mode  wave  of  the 
an^)htude 


|Cj(ii:,+0)|«|c,(i?.-0)|  (5.12) 

In  a  similar  way,  the  entropy  wave  of  the  an^litude  (R^  -  0)  generates  the  first  mode 

wave  of  the  amplitude 


|Cj(/?.+0)|«|cXK.-0)| 


(5.13) 
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Further  downstream  (see  Region  3  in  Figures  2  and  5)  the  first  mode  is  synchronized  with 
the  second  mode.  As  it  was  shown  in  Section  4,  the  first  mode  wave  excites  the  second 
mode  wave  of  the  ait^litude 


|c,(i^,+0)|«|c,(;?^i-0)|  (5.14) 

where  is  the  branch  point  located  in  Region  3.  Since  the  first  mode  attenuation  is 
relatively  weak  in  a  distance  between  points  and  ,  we  can  estunate  the  first-mode 
wave  an^htude  as 


lCi(i^i-0)|«|Ci(i?.+0)|  (5.15) 

Finally,  using  Eqs.  (5.6),  (5.7),  (5.11)  and  (5.12)  we  estimate  the  initial  amplitude  of  the 
normal  mode  waves  as 

\c,(R,,+0)i^cM,+0)  «|c,(^  -0)1  or  +0)\^\cM,  +0)  «|c,(7?.  -  0)|  (5.16) 

where  c„  and  are  anq)litudes  of  the  vorticity  and  entropy  wave  satisfying  the  synchronism 
condition  (5.4)  in  Region  2.  As  shown  in  Section  4,  one  of  these  normal  waves  gets 
essentially  unstable  downstream  fi'omthe  branch  point  .  This  is  the  Mack’s  second-mode 
instabihty.  Its  initial  anqtlitude  at  the  branch  point  can  be  roughly  estimated  using  the 
relations  (5.16).  These  estimates  show  that  the  vorticity  (entropy)  waves  of  a  firequency  a 
and  the  wavenumber  (a^),  which  satisfy  the  synchronism  condition  (5.5),  are 
“swallowed”  by  the  boundary  layer  in  the  synchronism  Region  2  and  generates  the  unstable 
normal  waves  of  the  same  fi'equency.  Sketch  of  this  swallowing  mechanism  is  shown  in  Fig. 
19.  The  unstable-wave  initial  amplitude  of  some  physical  value,  measured  at  the 

branch  point  i?  =  Re(i?ji(<a)),  is  approximately  equal  to  the  anq)litude  of  said  vorticity 
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(entropy)  wave,  ),  measured  near  the  synchronism  point  R,  (co)  -  0 .  In  wind  tunnel 

experiments,  vorticity  (entropy)  wave  amphtude  near  the  synchronism  point  R,,  (Qe), 

can  be  measured  using  a  hot-wire  anemometry  with  streamwise  and  vertical  probe 
separation  [42]  or  hot  films  placed  on  the  wall  near  the  sj^chronism  point  R .  Coupling 
between  vorticity/entropy  waves  and  the  Mack’s  second  mode  can  be  verified  by  correlation 
measurements. 

It  is  generally  believed  that  for  the  case  of  hypersonic  wind  tunnels,  the  major  source  of 
fi’eestream  disturbances  is  fi’om  the  radiation  of  sound  generated  by  the  turbulent  boundary 
layer  on  the  wind  tunnel  nozzle  walls.  However,  most  of  the  acoustic  energy  is  concentrated 
at  low  fi'equencies,  and  sound  amplitudes  decreases  with  increasing  frequency.  For  thin 
boundary  layers  at  high  unit  Reynolds  numbers,  the  second-mode  waves  may  have 
frequencies,  which  are  larger  than  the  fi'equencies  where  the  fireestream  sound  is  lost  in  the 
electronic  noise.  Therefore,  a  conventional  hypersonic  wind  tunnel  may  be  “quiet”  for  the 
second  mode  disturbances  [26],  In  this  case,  the  fi-eestream  entropy  and  vorticity 
disturbances  may  play  a  major  role  in  the  boxmdary-layer  transition.  Similar  situation  may 
occur  in  flight  conditions  due  to  the  presence  of  atmospheric  turbulence  and  the  absence  of 
fi'eestream  acoustic  disturbances.  The  transition  onset  point  can  be  correlated  with  the 
fireestream  vorticity  and  entropy  disturbances  using  the  PSE  method  [19]  with  the  initial 
amphtude  distribution  of  the  Mack’s  second  mode  being  estimated  as 

Here  and  are  amplitudes  of  the  vorticity  and  entropy  wave  satisfying  the 
synchronism  condition  (5. 10)  in  Region  2. 

Foregoing  analysis  addresses  the  entropy/vorticity  waves,  which  are  strictly  satisfy  the 
synchronism  conditions,  and  does  not  account  for  neighboring  waves.  Additional  exchange 
between  normal  modes  and  the  vorticity/entropy  waves,  which  are  shghtly  de-tuned  fi’om  the 
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synchronism  ,  may  affect  the  initial  amplitude  of  unstable  normal  waves.  However,  rough 
estimates  (5.17)  can  be  usefiil  for  the  transition  onset  prediction  due  to  the  following 
reasons.  According  to  the  amplitude  method  [21],  the  transition  onset  point  x,,.  is 
determined  from  the  equation 


(5.18) 


where  -a.  is  spatial  growth  rate,  x^  is  locus  of  the  instability  growth  onset,  is  initial 
amplitude,  and  is  critical  amplitude  related  to  the  beginning  of  nonhnear  processes. 
Since  the  transition  point  depends  weakly  (as  a  logarithm)  on  the  initial  anq)litude  value, 
even  rough  estimates  of  A^  can  provide  sufBcient  accuracy  for  the  transition  onset 
prediction.  Such  estimates  can  be  economically  made  using  the  coupling  (5.17)  and 
integrated  into  stability  codes. 


6.  Summary  discussion  and  future  effort 

Transition  process  is  an  initial  boundary  value  problem,  which  requires  distributions  of  the 
normal  mode  amplitudes  at  the  instability  growth  onset.  These  initial  data  are  determined 
from  receptivity  analysis  comprising  the  identification  of  external  disturbances,  the 
determination  of  the  most  sensitive  regions,  and  prediction  of  the  normal  mode  propagation 
from  these  regions  to  the  instability  growth  onset  (instability  pre-history).  Foregoing  analysis 
shows  that  all  these  elements  are  different  from  those  typical  for  subsonic  and  moderate 
supersonic  speeds.  Low  supersonic  stability  concepts  cannot  be  just  extended  to  hypersonic 
cases.  The  multiple-mode  method  rather  than  a  single-mode  approach  is  required  to  evaluate 
receptivity  and  instability  pre-history  of  hypersonic  boundary  layers.  These  phases  of  the 
transition  process  are  more  cottq)licated  than  in  moderate  supersonic  cases  due  to  the 
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presence  of  the  disturbance  spectrum  singularities.  As  shown  in  Figures  3a  and  6a,  these 
singularities  are  located  in: 

•  the  nose  region  1,  where  the  first  and  second  waves  are  synchronized  with  external 
acoustic  waves  of  the  phase  q)eed  =  \  +  \l  and  =  1  - 1  /  respectively. 

•  the  synchronism  region  2,  where  the  first  mode  is  synchronized  with  the  external 
entropy/vorticity  waves  of  the  phase  velocity  =  1 .  In  this  region,  the  first  mode 
penetrates  outside  the  boundary  layer.  Its  oscillating  “tail”  (see  Fig.  5)  could  be  identified 
by  measuring  correlations  between  hot-wire  signals  inside  and  outside  the  boundary 
layer. 

•  the  synchronism  region  3,  where  the  first  and  second  modes  are  synchronized  near  the 
branch  point  of  discrete  spectrum  In  this  region,  the  boundary-layer  disturbance  field 
has  abnormal  behavior,  which  could  be  detected  using  the  hot-wire  measurements. 

In  the  nose  region  1,  acoustic  waves  can  effectively  excite  the  first  and  second  modes  [32], 
[35].  This  receptivity  mechanism  is  quahtatively  different  fi’om  the  subsonic  case  reviewed 
by  Choudhari  [18].  In  conventional  hypersonic  wind  tunnels,  fi'eestream  acoustic  waves  are 
generated  by  sources  propagating  in  the  turbulent  boundary  layer  on  nozzle  walls.  These 
waves  have  the  phase  speed  c  <  1  - 1  /  ,  and  do  not  satisfy  the  synchronism  condition. 

However,  their  interaction  with  the  nose  tip  induces  a  local  source,  which  generates 
“dangerous”  disturbances  of  the  phase  speed  c  =  1  + 1  /  ,  as  schematically  shown  in  Fig. 

4.  For  a  planar  configuration  (flat  plate  or  hollow  cylinder),  this  acoustic  source  is 
concentrated  at  the  leading  edge.  For  the  conical  flow,  it  is  concentrated  near  the  cone  tip. 
Since  the  point  somce  is  weaker  than  the  line  source  and  downstream  attenuation  of 
acoustic  energy  fi-om  the  point  source  is  stronger  than  that  fi^om  the  line  source,  receptivity 
of  a  planar  configuration  to  acoustic  disturbances  should  be  higher  than  that  of  a  conical 
configuration.  These  arguments  are  consistent  with  experimental  observations  [33].  The 
fi^eestream  acoustic  scattered  by  sharp  leasing  edges  can  be  used  to  create  a  source 
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providing  disturbances  of  the  frequency  and  amplitude  rmder  control.  Schemes  of  such  a 
disturbance  generator  are  shown  iti  Figures  20a  and  20b.  When  an  external  acoustic  source 
(A)  is  placed  below  the  model  axis  (see  Fig.  20a),  then  it  induces  an  acoustic  source  of  the 
same  frequency  at  the  leading  edge.  Acoustic  disturbances  from  the  leading  edge  source 
penetrate  into  the  boundary  layer,  whereas  the  other  portion  of  disturbances  is  blocked  by 
the  body  and  do  not  affect  the  boundary  layer.  As  contrasted  to  natural  transition 
experiments,  this  scheme  allows  direct  measurements  of  the  hypersonic  leading  edge 
receptivity  and  downstream  evolution  of  the  boundary-layer  disturbances.  If  the  external 
acoustic  source  (B)  is  placed  above  the  model  as  shown  in  Fig.  20b,  then  the  boundary  layer 
is  radiated  by  acoustic  waves  from  both  Source  (B)  and  the  leading  edge  source.  Comparing 
these  two  cases  we  could  estimate  the  leading  edge  receptivity  versus  the  distributed 
receptivity. 

In  the  synchronism  region  2,  vorticity  and  entropy  waves  are  partially  swallowed  by  the 
bormdary  layer  and  effectively  generate  the  first  mode  as  schematically  shown  in  Fig.  19. 
This  receptivity  mechanism  can  conqpete  with  the  leading  edge  receptivity  to  freestream 
acoustic  disturbances  in  cases  wdien: 

•  acoustic  disturbances  are  small  compared  to  freestream  turbulence  and  tenq)erature 
spottiness  in  the  frequency  band  of  the  bormdary-layer  instability.  Low  level  of 
freestream  noise  is  typical  for  flight  conditions.  Moreover,  Stetson  and  Kimmel  [26] 
noted  that  for  the  wind  tunnel  experiments  at  high  Reynolds  numbers,  the  bormdary- 
layer  instability  may  have  frequencies,  which  are  larger  than  frequencies  where  the 
freestream  noise  is  lost  in  the  electronic  noise.  Therefore,  even  conventional  hypersonic 
wind  tunnels  may  be  “quiet”  in  the  instability  frequency  band. 

•  the  leading  edge  shape  corresponds  to  relatively  weak  interaction  between  external 
acoustic  field  and  the  boxmdary  layer.  The  sharp  cone  is  a  good  exan^le  of  this  case.  In 
contrast  to  the  conical  flow,  a  planar  configuration  (flat  plate  or  hollow  cylinder)  can 
create  a  strong  local  source  at  the  leading  edge  due  to  external  acoustic  scattering.  Since 
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the  point  source  on  the  cone  tip  is  weaker  than  the  line  source  on  a  plate  leading  edge 
and  downstream  attenuation  of  the  acoustic  energy  from  the  point  source  is  stronger 
than  that  from  the  line  source,  receptivity  of  a  planar  configuration  to  acoustic 
distiubances  should  be  higher  than  that  of  a  conical  configuration.  These  arguments  are 
consistent  with  experimental  observations  [33], 

Analysis  of  Section  5  provides  rough  estimates  for  coupling  between  the  external 
vorticity/entropy  disturbances  and  the  initial  an^litude  of  hypersonic  boundary-layer 
instability.  These  results  can  be  used  as  initial  data  for  the  PSE  analysis  downstream  from  the 
synchronism  region  2.  Since  the  transition  point  depends  weakly  (as  a  logarithm)  on  the 
initial  amplitude  of  the  boundary-layer  instability,  even  rough  estimates  of  the  receptivity 
mechanism  can  provide  sufficient  accuracy  for  the  transition  onset  prediction.  Such 
estimates  can  be  economically  made  using  the  coupling  relation  (5.17)  and  integrated  into 
stability  codes  such  as  the  PSE  code.  That  makes  feasible  an  amplitude  method  for  the 
transition  prediction  couplmg  receptivity  and  stability  problems. 

In  the  synchronism  region  3,  characteristics  of  the  first  mode  are  very  close  to  those  of  the 
second  mode.  The  discrete  spectrum  has  the  branch  point  located  in  the  conq)lex  R- 
plane  near  the  real  axis.  Near  the  branch  point  the  discrete  spectrum  sphts  into  two  branches, 
as  schematically  shown  in  Fig.  7.  Instability  of  the  Mack’s  second  mode  is  associated  with 
this  spHtting.  In  the  vicinity  of  the  branch  point,  the  normal  mode  decon^osition  is  not  vahd 
and  should  be  replaced  by  local  solution,  which  provides  coupling  between  input  and  output 
ampHtudes.  In  the  first  order  approximation,  the  inter-mode  exchange  rule  (4.21)  does  not 
depend  on  the  mean  flow  profiles  and  can  be  applied  for  broad  class  of  shear  layers  with 
similar  topology  of  the  discrete  spectrum. 

Calculations  of  the  branch  point  as  a  fimction  of  the  disturbance  frequency,  local  Mach 
number  and  wall  temperature  ratio  show  that  the  branch  pomt  can  intersect  the  real  axis  (see 
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Fig.  13b).  That  leads  to  different  topology  of  the  discrete  spectrum  splitting.  If  the  branch 
point  is  in  the  upper  half  of  the  complex  -plane,  it  is  bypassed  from  lower  side  as  the 
Reynolds  number  increases  along  the  real  axis.  Downstream  from  the  branch  point,  the  first 
mode  gets  essentially  unstable  and  the  second  mode  gets  more  stable.  This  case  is  shown  in 
Figures  6a  and  6b  and  corresponds  to  the  cooled  wall  cases  of  relatively  low  Mach  numbers, 
<  6 .  If  the  branch  point  is  in  the  lower  half  plane,  it  is  bypassed  from  the  upper  side. 
Downstream  from  the  branch  point,  the  second  mode  gets  unstable  and  the  first  mode  gets 
more  stable.  This  case  is  shown  in  Figures  3a  and  3b.  Infinitesimal  variation  of  the  basic 
parameter  (say  the  Mach  number  A/^)  near  the  value  at  which  the  branch  point  crosses  the 
real  axis  causes  jurt^  from  one  topology  to  another.  Moreover,  near  the  branch  point  the 
normal  wave  amplitude  tends  to  infinity  and  the  group  velocity  has  non-monotonic  behavior. 
These  features  along  with  the  spectrum  singularity  in  the  synchronism  region  2  can  cause 
significant  difficulties  in  the  PSE  integration  as  well  as  stability  calculations  using  the 
or  COSAL  codes. 

Summarizing  we  can  formulate  the  following  possible  scenario  of  the  transition  initial  phase: 

External  entropy/vorticity  disturbances  generate  the  first  mode  in  the  synchronism  region 

2.  This  mode  excites  the  second  mode  near  the  branch  point  of  the  synchronism  region 

3.  Receptivity  and  the  inter-mode  exchange  mechanisms  relevant  to  these  regions  provide 
initial  amplitude  for  the  Mack’s  second  mode  just  near  the  instability  growth  onset.  Further 
downstream,  the  Mack ’s  second  mode  amplifies  up  to  the  critical  amplitude  corresponding 
the  transition  onset. 

This  scenario  is  more  probable  for  the  conical  configurations  and  “quiet”  freestream 
conditions.  It  is  consistent  with  the  following  experimental  observations: 

•  Transition  on  a  sharp  cone  is  less  sensitive  to  the  cone  tip  conditions  (wall  temperature 
jiunp,  small  bluntness)  than  transition  on  a  flat  plate. 
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•  Hypersonic  cone  transition  data  [25]  do  not  fit  to  the  correlation  of  Pate  and  Scheuler 
based  on  assumption  that  external  acoustic  waves  trip  the  boundary  layer  (see  [44]). 

•  Since  the  Mack’s  second-mode  waves  are  excited  just  near  the  instability  onset,  they 
start  to  grow  with  the  initial  amplitude  to  be  relatively  slow  fimction  of  the  disturbance 
frequency.  These  waves  become  “tuned”  to  the  boxmdary  layer  thickness  and  tend  to  the 
optimal  tuning  near  the  envelope  curve.  That  is  consistent  with  experimental  data 
reported  in  [26]. 

This  study  is  only  the  first  step  in  understanding  of  hypersonic  boundary-layer  receptivity 
and  instabihty  pre-history.  However,  it  provides  good  starting  conditions  to  attack  the 
following  challenging  problems: 

•  Foregoing  analysis  should  be  extended  to  the  three-dimensional  disturbance  case.  Wave 
packets  rather  than  plane  waves  need  to  be  addressed.  That  will  help  to  evaluate 
receptivity  to  3-D  fi’eestream  disturbances  and  downstream  propagation  of  unstable 
wave  packets  typical  for  a  hypersonic  wind  tunnel  transition. 

•  Receptivity  to  vorticity/entropy  waves  needs  to  be  studied  in  detail.  Rough  estimates  of 
receptivity  mechanism  should  be  refined  from  systematic  analysis  of  the  unit  problems 
formulated  in  Section  5. 

•  The  flat-plate  leading  edge  receptivity  model  [32],  [35]  should  be  extended  to  the 
conical  configurations  and  the  blunt  leading  edge  cases. 

Investigation  of  these  problems  will  help  to  explain  hypersonic  stability  and  transition 
experiments  and  develop  reliable  transition  prediction  methods  required  for  hypersonic 
vehicle  design. 
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Fig.  1  Disturbance  spectrum  topology  for  supersonic  boundary  layer. 
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Fig.  3a  Disturbance  spectrum  for  the  conical  boxmdary  layer  at  F  =  10 
7^  =  7^ ;  the  eigenvalue  real  part  distributions  a  (R) . 


Fig.  3b  Disturbance  spectrum  for  the  conical  boundary  layer  at  F  =  10 
=  7^ ;  the  eigenvalue  imaginary  part  distributions  a. (R) . 


Fig.  4  Schematic  pattern  of  interaction  between  an  external  acoustic  wave  and  the  leading 
edge. 


Fig.  5  The  first  mode  eigenfimction  for  the  conical  boundary  layer  at  F  =  10^ ;  =  6.8 , 

T  =  T 

-'h-  -'orf  • 
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Fig.  6a  Disturbance  spectrum  for  the  conical  boundary  layer  at  F'  =  10^ ;  =  5.5 , 

^  ^  =  0.1 ;  the  eigenvalue  real  part  distributions  a ^{R) . 


Fig.  6b  Disturbance  spectrum  for  the  conical  boundary  layer  at  F’  =  10"^ ;  =  5.5 , 

=  0.1 ;  the  eigenvalue  imaginary  part  distributions 
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Fig.  8  Stokes  lines  Ij  and  canonical  domains  Dj  in  the  branch  point  vicinity 


unstable  unstable 


stable  stable 


(a) 


(b) 


Fig.  9  The  inter-mode  exchange  diagrams;  (a)  -  stable  mode  transformation;  (b)  -  unstable 
mode  transformation. 


61 


•  •  -  Dallmann,  conical  metric 
A  A  -  Stetson 

-a,  lO’ 


Fig.  11  Location  of  the  branch  point  compared  to  experimental  data  of  Stetson  et  al. 
[25]  and  linear  stabihty  calculations  of  Simen  and  Dallmann  [28];  A/^=6.8, 
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Fig.  12  Location  of  the  branch  frequency  compared  to  e)q)ernnental  data  of  Stetson  et  al. 
[25]  and  linear  stability  calculations  of  Simen  and  Dallmann  [28];  M^=6.8, 

T^  =  T^,R=m. 
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Re(Rb,)10'^ 


Fig.  13a  Real  part  of  the  branch  point  as  a  function  of  the  local  Mach  number  at 

various  wall  tenq)erature  ratios  IT^',  F  =  10^ . 


Tm(Rbi)10-^ 


Fig.  13b  Imaginary  part  of  the  branch  point  as  a  function  of  the  local  Mach  number 

at  various  wall  temperature  ratios  T^IT^  \  F  =  10^ . 
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Fig.  14a  Disturbance  spectrum  for  the  conical  boundaiy  layer  at  F  =  10"^ ;  =  5.95 , 

/  7^  =  0.1 ;  the  eigenvalue  real  part  distributions  . 


Fig.  14b  Disturbance  spectrum  for  the  conical  boimdary  layer  at  F’  =  10^ ;  =  5.95 , 

7^  /  7^  =  0.1 ;  the  eigenvalue  imaginary  part  distributions  a,.  (7?) . 
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Fig.  15a  Disturbance  spectrum  cc^{R)  in  the  vicinity  of  synchronism  point  R,  w  1515 ; 
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Fig.  16  Spectrum  topology  in  the  complex  a  -plane  near  the  ^uchronism  point  R, . 


Fig.  17  First  mode  eigenfunctions  of  the  x-con^onent  velocity  u(tj)  for  various  Reynolds 

numbers;!  -  /?=  1369;2-  7?  =  1438;  3-  7?  =  1507;  4-  7?  =1542;  5-  7?  =  1611;6- 
7?  =1681. 
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disturbance  swallowing 


R,(F,)  R,(Fp 


Fig.  19  Schematic  picture  of  the  freestream  vorticity/entropy  disturbance  swaflowing  in 
synchronism  regions  near  the  points  R, (F^ )  and  ;  frequency  F^  >  F^. 
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Acoustic  field  from  the  leading  edge  source 


Fig.  20a  Scheme  of  experiment  with  the  acoustic  source  under  control;  external  source 
(A)  is  below  the  model. 


Acoustic  field  from  the  leading  edge  source 


Fig.  20b  Scheme  of  experiment  with  the  acoustic  source  under  control;  external  source 
(B)  is  above  the  model. 


71 


